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ABSTRACT 

Context. Near-Earth asteroid 162173 (1999 JU3) is a potential flyby and rendezvous target for interplanetary missions because of its 
easy to reach orbit. The physical and thermal properties of the asteroid are relevant for establishing the scientific mission goals and 
also important in the context of near-Earth object studies in general. 

Aims. Our goal was to derive key physical parameters such as shape, spin-vector, size, geometric albedo, and surface properties of 
162173 (1999 JU3). 

Methods. With three sets of published thermal observations (ground-based N-band, Akari IRC, Spitzer IRS), we applied a thermo- 
physical model to derive the radiometric properties of the asteroid. The calculations were performed for the full range of possible 
shape and spin-vector solutions derived from the available sample of visual lightcurve observations. 

Results. The near-Earth asteroid 162173 (1999 JU3) has an effective diameter of 0.87 ±0.03 km and a geometric albedo of 
0.070 ±0.006. The^ 2 -test reveals a strong preference for a retrograde sense of rotation with a spin-axis orientation of /t ed = 73°, 
y8c C i = -62° and P s id = 7.63±0.01 h. The most likely thermal inertia ranges between 200 and 600 JirT 2 s~ 5 K~', about a factor 
of 2 lower than the value for 25143 Itokawa. This indicates that the surface lies somewhere between a thick-dust regolith and a 
rock/boulder/cm-sized, gravel-dominated surface like that of 25143 Itokawa. Our analysis represents the first time that shape and 
spin-vector information has been derived from a combined data set of visual lightcurves (reflected light) and mid-infrared photometry 
and spectroscopy (thermal emission). 

Key words. Minor planets, asteroids: individual - Radiation mechanisms: Thermal - Techniques: photometric - Infrared: planetary 
systems - 162173 (1999 JU3) 



1. Introduction 



Asteroid 162173 (1999 JU3) is currently among the potential 
targets of future interplanetary exploration missions. The target 
is relatively easy to reach with state-of-the-art mission capabil- 
ities, and it offers high scientific potential (Binzel et al. 2004). 
This near-Earth asteroid belongs to the C-class objects, which 
are believed to represent primitive, volatile-rich remnants of the 
early solar system. Various aspects of this small body have been 
covered in some detail in the recent works by Hasegawa et al. 
d20Q8l i and by Campins et al. (120091 . 



Hasegawa et al. (20Q8J use a spherical shape model for their 
radiometric analysis, and alternatively use an ellipsoidal shape 
model, but without knowing the true spin-vector orientation. 
The results (radiometric diameter of 0.92 + 0.12 km, visual ge- 
ometric albedo of 0.063^^5) which indicate a thermal inertia 
larger than 500 J m~ 2 s~ a5 K -1 ), were based on a set of photomet- 
ric Subaru and Akari observations and connected to simplified 
shape and spin-axis assumptions. Campins et al. (2009) have ob- 
tained a single-epoch Spitzer infrared spectrograph (IRS) spec- 
trum. They used a spherical shape model, and for the spin-pole 
orientation they used the extreme case of an equatorial retro- 
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grade geometry and a prograde solution published by Abe et al. 
(2008). Their analysis, based on the single IRS-spectrum and ig- 
noring the data sets published by Hasegawa et al. (2008), yielded 
a value for the thermal inertia of 700 + 200 J m -2 s~°- 5 K -1 , a 
diameter estimate of 0.90 + 0.14km, and geometric albedo of 
0.07 + 0.01. 

Despite the simplifications in shape and spin-vector prop- 
erties, both sets of published radiometric diameter and albedo 
values agree within the given uncertainties. Both teams also 
favoured a relatively high thermal inertia (close to that of 
25143 Itokawa). The unknowns of the spin- vector orientation 
cause a large uncertainty in the thermal properties. 

Here we re-analyse all available lightcurve observations to 
derive (on the basis of standard y 1 lightcurve inversion tech- 
niques) all matching spin-vector and shape solutions (Sect. [2]). 
The full possible range for shape, spin-axis orientation, and rota- 
tion period was then used as input for a thermophysical ;^ 2 anal- 
ysis of all available thermal observations (Sect. [3]) with the goal 
of deriving radiometric sizes, albedos and thermal inertias. At 
the same time, we determined the most likely shape-solution, 
rotation period and spin-axis orientation (Sect. [4]). 



2. Possible shape and spin-vector solutions 

A detailed list of the available photometric observations is pre- 
sented in TableQ] There are about 40 dedicated visual lightcurve 
data sets spread over more than 270 days. About half of the 
lightcurves were calibrated; some of the data sets are very noisy. 
Based on these data, Abe et al. ((2008 ) found a rotation period 
of 7.6272 + 0.0072 h and a spin orientation of /l ec i = 331.0°, 
jSeci = +20.0°, indicating a prograde rotation. A more recent 
analysis by the same authors (priv. communication) resulted in 
a rotation model with slightly different values: /l ec i = 327.3°, 
jSed = +34.7°, P sid = 7.6273922 hours. Both solutions were 
derived using the epoch and amplitude methods described by 
Magnusson (1986). These methods are reliable for irregularly 
shaped bodies and for sufficient lightcurve data covering vari- 
ous aspect angles. 1999 JU3 has a comparatively spherical shape 
and the available lightcurve data were from restricted directions. 
Nevertheless, both solutions gave a reasonable match to the ob- 
served lightcurves, but it turned out that these solutions are not 
unique and other parameter sets with different rotation periods, 
spin-axis orientations and shape models could not be excluded 
(see Fig.Q}. 



Table 1. Observation circumstances for the lightcurve measure- 
ments (see also Table. 1 in Abe et al. 2008). 



Mon/Day (2007) 



Telescope 



Observer 



07/8, 09/4 

07/19-23, 12/3,4,6-8, 
02/26-28, 04/2,4,5 

08/5,15,09/6,11,13,15, 
10/6,18, 11/13,15 

08/9-10,17,20, 09/6,10 

09/4,5,7,8,10,12,14,15 
11/7-9,11,13, 02/5-8, 
04/14,15 

09/11-14 



2.2 m/Mauna Kea 
l.Om/Lulin 

l.Om/Ishigaki 

1.0 m/Bisei 

1.05m/Kiso 

1.55 m/Steward 



T. Kasuga 

M. Abe, K. Kawakami, 
D. Kinoshita 

D. Kuroda, S. Nagayama, 
K. Yanagisawa 

S. Urakawa, 
S. Okumura 

M. Abe, K. Kawakami, 
Y. Sarugaku, Y. Takagi, 
S. Miyasaka 

P. R. Weissman, 
Y.-J. Choi, S. Larson 



Since these values are crucial input parameters for our ther- 
mophysical model analysis, we repeated the search for possible 
shape and spin-vector solutions using the lightcurve inversion 
method developed by Kaasalainen and Torppa (2001 ). Our goal 
was to derive a set of the most likely convex shape models that 
would fit all available lightcurves. Because of the poor quality 
of some of the data, we were only able to determine the range 
of the sidereal rotation period to 7.6204-7.65 lOhours. Different 
shape models with pole orientations covering almost the entire 
celestial sphere (without any preference for pro- or retrograde 
solutions) fit the data equally well. By scanning the period-pole 
parameter space, we derived 77 shape models that were physi- 
cally acceptable (rotating around the shortest axis) and for which 
fhe^f 2 of the fit was no more than 10% higher than the best-fit 
X 2 . These models corresponded to local minima in the parameter 
space. 

In addition, we included both of the original Abe et al. (120081 ) 
solutions and added another five with the pole fixed to the two 
(/led, /8eci)-pairs mentioned above and rotation periods in the 
range given by Abe et al. (2008). Only two of these models 
pass the^ 2 + 10% limit on the basis of the visual lightcurves, 
the other five models have higher^- 2 -values. Six of these shape 
models are elongated along the 'z' axis and therefore unphysical. 

For all 84 models we performed a thermophysical model 
analysis for a wide range of possible parameters (see Table|3). 



2007/7/23.7 2007/9/13.6 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

2008/2/28.5 2008/4/5.6 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

Phase of rotation Phase of rotation 



Fig. 1. Match between observed lightcurves and shape/spin-axis 
model solution "7_1". 



3. Thermophysical Model Analysis 

3.1. Model and input parameters 

The mid-IR photometric data were already described in 
Hasegawa et al. (12008ft . The data set includes 15 N-band Subaru 
observations and two dedicated Akari observations at 15 and 
24/im. We binned the single-epoch Spitzer IRS data (Campins 
et al. 120091 into 20 wavelength points (4 for band SL2, 7 for 
SL1, 4 for LL2, 5 for LL1; see Fig. |5] bottom). The 20 wave- 
length points were chosen to give about equal weight to the two 
published data samples in terms of number of observations (17 in 
Hasegawa et al. 2008 and 20 for the Campins et al. |2009l sample). 
In this way the derived object properties are better connected to 
the entire data set, and they do not just match the measurements 
of one observing epoch. Each observation set also has a mixture 
of lower and higher quality data: The ground-based Subaru data 
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are of lower quality than the Akari data, while the quality of the 
IRS spectrum changes with wavelength. This again ensures that 
the final solutions are not biased towards a single measurement. 
All observations are listed in Table|2] 

All 84 possible spin-vector and shape solutions from Sect. [2] 
have been used in combination with these thermal data. 

For our analysis we are using a thermophysical model 
(TPM) de scribe d by Lagerros ([19961 [19971 [19981 and Miiller & 
Lagerros ( 1998 ). This TPM works with true illumination and ob- 
serving geometries, accepts irregular shape models and arbitrary 
spin-vector solutions, works with roughness controlled by the 
r.m.s. of the surface slopes, considers heat-conduction into the 
surface as well as multiple scattering of both the solar and the 
thermally emitted radiation. The model has been tested and vali- 
dated thoroughly for NEAs (e.g., Miiller et al. 120051 ) and MBAs 
(e.g., Miiller & Lagerros 2002). The TPM input parameters and 
applied variations are listed in Table. [3] 



Table 3. Summary of general TPM input parameters and applied 
ranges. 



Param. 


Value/Range 


Remarks 


r 


0...2500 


Jm~ 2 s -0 ' 5 K" 1 , thermal inertia 


P 


0.4...0.9 


r.m.s. of the surface slopes 


f 


0.4...0.9 


surface fraction covered by craters 


e 


0.9 


emissivity 


i/v-mag. 


18.82 ±0.02 mag 


Abeetal. (2008) 


G-slope 


-0.1 10 ±0.007 


Abe etal. (2008) 


shape 


84 models 


see Sect. [2] 


spin-axis 


84 solutions 


see Sect. [2] 


P s ,d [h] 


7.6205. ..7.65 10 


see Sect. [2] 



3.2. Solving for effective diameter, geometric albedo and 
thermal inertia 

Campins et al. (2009) and Mueller (120071 ) used a^ 2 or reduced 
^ 2 -test to find solutions for the thermal inertia F. Here we fol- 
low a modified approach to find the most robust solutions with 
respect to thermal inertia and allowing for the full range in ef- 
fective diameter and geometric albedo at the same time. The fol- 
lowing procedure was executed for all 84 possible shape and 
spin- vector solutions separately: 

(1) For each value of T in a wide range (see Table [3} we calcu- 
late the radiometric diameter and albedo solution via the TPM 
for each individually observed thermal flux (37 individual diam- 
eter and albedo solutions). Diameter and albedo are linked by the 
absolute magnitude Hy which was kept constant (the rotational 
amplitude is only about 0.1 mag). 

(2) We calculate the weighted mean radiometric diameter and 
albedo solution for each given F (x — ^'.""-i , with diame- 

ter/albedo errors cr, connected to the observational errors). 

(3) For each individual observation we predict TPM fluxes based 
on the given F and the corresponding weighted mean radiomet- 
ric diameter and albedo from step (2). 

(4) The most robust solutions occur when the observations and 
the TPM predictions agree best (taking the uncertainties of the 
measurements into account in a weighted mean sense, see step 
2). This can be expressed as jj'L^dobsi - mod^/cri) 2 , a modi- 
fied reduced^- 2 method. The most likely thermal inertia is found 



at the smallest x 2 values; the connected effective diameter and 
geometric albedo values are the ones calculated at step (2). 

In a first round of executing this procedure we kept the 
surface roughness constant at values which were specified as 
"default values" for large, regolith-covered asteroids (Miiller et 
al. [1999). The corresponding roughness parameter values are 
p-0.7, the r.m.s. of the surface slopes, and f-0.6, the fraction 
of surface covered by craters. The results are shown in Fig. [2] 




10 100 1000 

Thermal inertia [J m 2 s 05 K '] 



Fig. 2. TPM x 1 -optimization process to find robust solutions for 
diameter, albedo and thermal inertia simultaneously. Each line 
represents the reduced x 1 values for an individual shape/spin- 
vector solution as a function of thermal inertia. The surface 
roughness was kept constant using the "default values" of p=0.7 
and /=0.6 as proposed by Miiller et al. (119991) . Model 7_1 solu- 
tions marked with squares on the solid line correspond to the two 
cases with thermal inertias of 20 and 1000 J m -2 s~ K l shown 
in Fig. [3] (left and middle). 



The best shape and spin-vector solutions (with lowest val- 
ues for the reduced x 2 and clear minima in Fig. [2| were then 
the starting point for further tests: (i) Are these solutions robust 
against sub-sets of the thermal data? (ii) How does the surface 
roughness influence the solutions? (iii) Do the solutions explain 
the thermal behaviour over the observed phase angle range (from 
-20° to -55°)? (iv) Is the TPM match of equal quality at all ob- 
served wavelengths? (v) Are there large discrepancies at certain 
rotational phases? 

4. Results and Discussion 

4.1. Solution for shape and spin-vector 

The shape and spin-vector solutions which produce the lowest 
^-values in Fig.[2]are listed in Table|4] The model IDs represent 
a full shape-model, each with more than 2000 surface elements 
and more than 1000 vertices. The Julian date at zero rotational 
phase yo is m all cases 7o=2454289.0. 

All three solutions are retrograde solutions and, in fact, the 
eight best solutions in Fig. [2] are retrograde solutions. The best 
prograde solutions in the^ 2 picture are models with ID 5_2 and 
ID 13_5, both have ^-minima at 2.3 (almost a factor of 4 higher 
than the best retrograde solution 7_1) and would require an ex- 
tremely high thermal inertia (> 1000 J irT 2 s~ - 5 K _1 ) to match 
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Table 2. Summary of the avaible thermal observations. R/, is the helio-centric distance and A the distance between object and 
telescope. All observations were taken at positive phase angles a (Sun-object-telescope), i.e., leading the Sun, with a cold terminator 
for a retrograde rotating body. 



Year/Mon/Day Wavelength range [pm] Rj, [AU] A [AU] a [°] Telescope Reference 

2007/05/16 15.0,24.0 1.414 0.992 +45.6 Akari Hasegawa et al. J2008l > 

2007/08/28 8.8(3x), 9.7(lx), 10.5(lx), 11.7(7x), 12.4(3x) 1.287 0.306 +22.3 Subaru Hasegawa et al . f2008l > 

2008/05/02 5.2-8.5,7.4-14.2,14.0-21.5,19.5-38.0 1.202 0.416 +52.6° Spitzer Campins et al. ( |2009] > 



Table 4. The shape and spin-vector solutions which produce the 
lowest ^-values in Fig. [2] 



model-ID 


^ecl [°] 


ftol [°] 


P sid [h] 


^ 2 -min 


7.1" 


73.1 


-62.3 


7.6323 


0.60 


8.2" 


69.6 


-56.7 


7.6325 


0.68 


14.8 


77.1 


-30.9 


7.6510 


1.46 



Notes. (a) Models 7_1 and 8.2 are in the same local minimum in the 
parameter space for the lightcurve fits. 



the observations. Our most likely solution can be summarized 
as: A ecl = 73°, /3 e ci = -62°, P sid = 7.63 + 0.01 hours. Fig. [3] (left, 
middle) shows the model-ID 7.1, as seen from Spitzer during 
the IRS-observations and for the two thermal inertias marked 
in Fig. |2] with boxes. The match between model "7.1" with ob- 
served lightcurves is shown in Fig.Q] 

We also analysed the thermal data set against the Abe et al. 
(2008) spin-pole solutions discussed above. The corresponding 
^-minima are more than a factor five worse than our two best 
models above (see dashed lines in Fig. [2]). Some of these so- 
lutions are compatiple with the inertia range given by Campins 
et al. (120091 ) and even produced an excellent match to the IRS- 
spectrum. Nevertheless, these solutions can be excluded with 
very high confidence: (i) the match to the rest of the data set 
(Akari and ground-based data) is very poor (reflected in the high 
^-values); (ii) the corresponding shape models are unphysi- 
cal with elongations along the spin axis. Such rotational states 
would not be stable. The best of these models (in terms of x 2 - 
minima) is shown in Fig. [3] on the right side (with the rotation 
along the z-axis with the object's largest extension). 

4.2. Thermal inertia and surface roughness 

The thermal inertia is clearly a key parameter when modelling 
the mid-IR data for NEAs; it strongly influences the shape of 
the spectral energy distribution (SED). This can be seen in the 
Spitzer IRS spectrum, especially in the Wien-part of the spec- 
trum. But the thermal inertia also drives the thermal behaviour 
as a function of phase angle (e.g., Miiller 2002), the thermal 
phase curve. The temperature of the unilluminated fraction of 
the surface changes strongly with thermal inertia. In Fig. [3] the 
unilluminated fraction has just rotated out of the solar insola- 
tion. In one case (left) we assumed a low thermal inertia of 
~20 J rrT 2 s~ - 5 K _I , which was found to be typical for large 
MBAs (Miiller & Lagerros l2002l l. while in the second case (mid- 
dle) the thermal inertia is 1000 J m -2 s~ 5 K _I , close to the value 
found by Miiller et al. d2003l) for 25143 Itokawa. 



The importance of the thermal inertia in the modelling is also 
visible in the ^-solutions (Figs. [2] & @): the ^-values change 
significantly when going through the whole grid of physically 
meaningful thermal inertias. 



10 



5 



3 best shape/SV solutions 




10 



1000 



Fig. 4. TPM ^-optimization process for the shape and spin-axis 
solutions with the lowest ^ 2 -values (models 7_1, 8_2, 14.8). The 
solid lines are the "default roughness" values derived for MBAs. 



But Fig. |4] demonstrates that thermal inertia and surface 
roughness are not easy to disentangle, at least on the basis of 
this data set. Both parameters influence the short-wavelength 
SED-part where the hottest surface temperatures dominate the 
SED-shape. One way of solving for both parameters would be 
to obtain data with a larger phase angle coverage and larger 
wavelength coverage. The roughness plays a much bigger role 
at small phase angles (beaming effect), but thermal data close to 
opposition are not available for 162173 (1999 JU3). The ther- 
mal inertia is more important at large phase angles and at longer 
wavelengths where the disk-averaged temperature dominates the 
SED shape. In general, the larger the coverage in phase angle and 
wavelength, the more accurate is the determination of diameter, 
albedo, thermal inertia and roughness. 

Fig-0also shows that similarly low^ 2 -values can be reached 
independent of the surface roughness. This demonstrates that 
roughness effects are still important for the interpretation of 
thermal data (via the thermal phase curves), even at these rel- 
atively large phase angles (Miiller 2002). But the data set does 
not allow us to constrain the surface roughness which broadens 
the possible range of thermal inertias as can be seen in Fig. [4] 
Based on our three best shape/spin- vector solutions and con- 
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I 

-2( 

-1; 

I 



Fig. 3. Left and middle: TPM implementation of the shape model 7_1 at the time of the Spitzer IRS observations on 2008-May-02 
02:01 UT, as seen from Spitzer in asteroid-coordinates, i.e., the z-axis goes along the rotation axis. Left: low thermal inertia. Middle: 
high thermal inertia. Both solutions are marked with squares in Fig. [2] Right: Shape model fixed on the Abe et al. (2008 ) spin-vector 
and tuned to match the Campins et al. (2009) findings. The rotation along the largest object extension (z-axis) is unphysical. The 
surface temperatures are given in Kelvin. 



sidering the uncertainties due to roughness, we conclude that 
the thermal inertia is in the range 200-600 J rrT 2 s _0i K . This 
range is in agreement with the lower limit of 150 J rrT 2 s~ a5 K _1 
given by Campins et al. (2009 ) but lower than their best fit value 
of 700 + 200 J irT 2 s~ 05 K _1 . Our value is about a factor of two 
lower than the one found for 25 143 Itokawa (Miiller et al. 2005 ). 
We expect that the surface of 162173 (1999 JU3) is therefore dif- 
ferent in the sense that there might be fewer rocks and boulders 
and that the surface includes millimetre sized particles (as op- 
posed to the cm-sized gravel on Itokawa). We also find a rigorous 
lower limit to the thermal inertia of about 100 J irT 2 s~ - 5 K _1 , 
similar to Campins et al. (|2p09). This limit is not compatible 
with a thick dusty regolith covering the entire surface which 
would result in a very low thermal conductivity and thermal in- 
ertia values well below 1 00 J rrT 2 s~ K _1 , which would in prin- 
ciple be possible for an object of that size and the relativly slow 
rotation rate. For comparison, the Moon's thick regolith gives a 
value below 40JirT 2 s" 5 Kr 1 (Keihm [19841 calculated for T = 
300 K). 



4.3. Radiometric diameter and albedo solution 

Our best fit to all observations, as represented in Fig. [2] by 
the solid line, resulted in a radiometric effective diameter 
of 0.87 + 0.02 km and 0.070 + 0.003 for the visual geometric 
albedo. Both values are within the error bars of the solutions 
found by Hasegawa et al. d2008 ) and by Campins et al. ([2009), 
but now with much smaller errors. The uncertainties are based 
on the best x 1 -values for model-IDs "7_1" and "8_2" and the full 
variation in roughness (as shown in Fig. [4]). The radiometric ef- 
fective diameter is connected to the most likely shape model and 
spin-vector solution from above and corresponds to the size of a 
spherical object with equal volume. Using the determined pos- 
sible thermal inertia range of 200-600 J m" 2 s -0 ' 5 K _I and uncer- 
tainties in roughness and Hy -magnitude lead to a possible value 
for the effective diameter of 0.87 ±0.03 km. The derived geo- 
metric albedo provides the best solution between all thermal ob- 
servations over the full wavelength and phase angle range and 
the absolute Hy -magnitude. The Z/v-magnitude was given with 
only 0.02 mag error. Using a more realistic Hy-enor of +0. 1 mag 
leads to a final geometric albedo value of 0.070 + 0.006. The 
small uncertainties reflect the importance of multi-epoch, multi- 



wavelength and large phase angle coverage for thermophysical 
studies of small bodies. In a similar study by Miiller et al. (2005 ), 
also based on a large thermal data set and a shape model from 
lightcurve inversion techniques, the derived effective diameter 
agreed within 2% of the true, in-situ diameter (Miiller et al. 
2005; Fujiwara et al. 2006). The quoted uncertainties above are 
formal errors from the x 2 optimization, including the possible 
range in thermal inertia, roughness and Hy. 

The remaining discrepancy between the Spitzer and the 
Akari-fiux at 15 /.im (seen in Fig.[5]bottom) might be related to an 
additional error introduced by the flux scaling done by Campins 
et al. (120091 ) to match the short-wavelength part of the spectrum 
(<14/mi) to the long-wavelength part of the spectrum (>14jum). 
The mismatch is caused mainly by the placement of the object 
within the IRS slit and reflected in the specified 10% systematic 
absolute calibration uncertainty. 

Although the Fig. [4] solutions and the model match in Fig. [5] 
look very convincing, there are some uncertainties remaining: 
The spin-vector and shape solution from lightcurve inversion 
techniques is not very robust; more lightcurve observations are 
needed to improve the quality. The^ 2 -test works best if the ther- 
mal observations cover a wide range of wavelengths, phase an- 
gles (before and after opposition) and rotational phases. But all 
thermal observations have been taken at pre-opposition (posi- 
tive) phase angles (leading the Sun) where the unilluminated 
part of the surface visible to the observer is either warm (pro- 
grade rotation) or cold (retrograde rotation); observations after 
opposition are not yet available. A combined before and after 
opposition data set would constrain the sense of rotation and the 
thermal properties much better, observations close to opposition 
would determine the surface roughness a bit better, hence con- 
strain the thermal inertia further. 

It is also important to note here that higher thermal inertias 
(>600 J irT 2 s~° 5 K _1 , Fig. |3]middle) would make a slightly bet- 
ter fit to the Spitzer IRS spectrum (improvement mainly at the 
shortest wavelengths in Fig. [5] bottom), but would cause a sig- 
nificant dependency in the diameter and albedo solutions with 
phase angle. The measurements taken at around 20° phase an- 
gle (Subaru) would then have fluxes that are about 40% higher 
than the corresponding model predictions (i.e., values >1.4 in 
Fig. [5] top). The Akari fluxes would still be ~ 20% higher than 
the model predictions. Our observational data set covering about 
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Fig. 5. All observations divided by the corresponding TPM pre- 
dictions based on our optimized radiometric solution. Top: as a 
function of phase angle. Bottom: as a function of wavelength. 
The full set of Spitzer IRS data are shown with little circles 
(Campins et al. 120091 1, the triangles are the re-binned data, the 
Akari IRC data are represented by squares; the Subaru-COMICS 
observations by diamond symbols. 



lent match of the flux changes with phase angle (the phase angle 
range covered here is from ~20° to ~55°). 

The example of 162173 (1999 JU3) shows that a combina- 
tion of visual lightcurves (reflected sunlight) and mid-/far-IR 
photometry or photo-spectroscopy (thermal emission) can im- 
prove the quality of shape and spin-vector solutions significantly. 
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30° in phase angle constrains the possible thermal inertia to val- 
ues below 700 J nr 2 s" - 5 KT 1 . 



5. Conclusions 

The radiometric analysis provides the following results: (i) a 
strong indication of a retrograde sense of rotation; (ii) a spin- 
vector with A ec i = 73°, /3 ec i = -62°, P sid = 7.63 ± 0.01 h and y = 
at To = 2454289.0, derived for the first time based on a com- 
bined analysis of visual lightcurve data and thermal observa- 
tions; (iii) a shape model (here labelled with 7_1) as shown in 
Fig.[3](left & middle); (iv) a thermal inertia in the range 200 to 
600JtrT 2 s~ - 5 K _1 ; (v) a radiometric effective diameter (of an 
equal volume sphere) of D e jj = 0.87 + 0.03 km; (vi) a radiomet- 
ric geometric albedo of py = 0.070 + 0.006; (vii) a lower thermal 
inertia than for Itokawa, suggesting the presence of smaller par- 
ticles, < cm-sized, in the regolith, though likely not fine dust; 
(viii) very good agreement in the radiometric solutions between 
the Spitzer, the Akari and the Subaru observations; (ix) an excel- 
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